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Beschreibuno 

Verfahren und Gerat zur adaptiven Merkmalsanderung 
bei eindimensionaien Signalen 

10 

Die Erflndung betrifft ein Signal verarbeitungsverfahren gemafl dem OberbegrifF des Patentanspruchs 
1 sowie ein Gerat zur Realisierung dieses Verfahrens, welches gemafl einem der Anspriiche 14 bis 19 
ausgebildet ist. 

In der Technik miissen informarionstragende zeitiiche Verlaufe einzelner physikalischer GroBen (z.B. 

15 Schalldmckschwankungen bei einem alaistischen Signal) haufig nicht nur verstarkt, sondern auch in 
ihren Eigenschaften verandert werden. Solche Verlaufe werden als eindimensionale Signale bezeich- 
net. Die angestrebte Veranderung der Eigenschaften solcher Signale wird hier verallgemeinemd als 
Merkmalsanderung bezeichnet. Wegen der damit verbundenen praktischen Vorteile erfolgt diese 
Merkmalsanderung heute in der Regel mit elektronischen Mitteln. Dies erfordert, daD das in seinen 

20 Merkmalen zu verandemde Signal eine elektrische Spannungsschwankung darstellt. Es wird durch 
einen geeigneten Wandier in eine solche Spannungsschwankung umgewandelt (z.B. durch ein Mikro- 
fon bei einem akustischen Signal). Nach der Merkmalsanderung ist das Signal haufig wieder in die 
ursprQngliche physikalische Darstellungsfonn zuruckzuwandeln, wofur ein geeigneter Wandier erfor- 
derlich ist (z.B. Lautsprecher bei einem akustischen Signal). Haufig werden fur die Merkmals- 

25 anderung digitale Verfahren verwendet, weii diese eine Reihe wohlbekannter Vorteile aufweisen. 

Oft wird eine dem Prinzip nach bekannte und als frequenzseiektive Filterung bezeichnete Merkmals- 
anderung angewandt. Hierbei werden gewisse Teiie des Signalspektrums hervorgehoben, andere abge- 
schwacht. Eine weitere vorbekannte Gruppe von Merkmalsandeningen wird als Optimalfilterung oder 
st6rungsreduzierende Filterung bezeichnet. Hierbei geht es darurn, Anteile des Signals, welche von 

30 StGrquellen oder unerwOnschten Signalpfaden (z.B. akustischen Rtickkopplungspfaden) herrtlhren, zu 
erkennen und deranig vom Signal abzutrennen, daQ das gewiinschte Nutzsignal moglichst unverfalscht 
zurilckgewonnen wird. Eine weitere bekannte Gruppe von Merkmalsandeningen wird als Kompres- 
sion bezeichnet. Hierbei sind die Signalamplituden in Abhangigkeit von ihrer mornentanen Intensitat 
mehr oder weniger abzuschwiichen. Beispiele sind die Ubertragung eines Signals hoher Dynamik uber 

2; einen Run'dfunkkanal mit begrenztem Modulationshub oder die Anpassung eines Sprachsignais an den 
verminderten Dynamikbereich eines geschadigten Gehors mittels eines geeigneten Horgerats, wie z.B. 
im US-Patent 3894195 (K.D. Kryter) beschrieben. Eine weitere Form von Merkmalsandeningen ist 
die als Preemphasis bekannte voriibergehende Anhebung bestimmter Signalanteile, mit dem Ziel, das 
Signal robuster gegen den EinfluQ von Storungseintragungen in einem Ubertragungskanal zu machen. 
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um N Abtastwerte verschobenen Wiederholungen des zu langen Faltungsergebnisses (sog. zirkulare 
Falcung), wie ebenfalls in dem o.g. Buch von Oppenheim und Schafer erlautert wird. Diese verscho- 
ben uberlagerten Signal werte werden als Fehler im Ausgangssignal horbar (sog. 'Time domain 
aliasing"). Um diese Fehler wenigstens zu verkieinern wurde bisher vorgeschlagen, die Lange von 
5 FFT und IFFT wesentlicb groBer zu wahlen als die Eingangsstgnal-Blocklange M, um moglichst 
wenige Uberlagerungen innerhalb der Grundperiode der Lange N zu erhalten. Dadurch wird der 
Rechenaufwand drastisch erhoht, ohne daB eine sichere Fehlerbegrenzung moglich wird. Auch die 
Anwendung eines weiteren Fensters auf die von der IFFT gelieferten Ausgangsdatenblocke, die 
bereits versucht wurde, erlaubt keine sichere Fehlerbegrenzung, fuhrt andererseics aber zu weiteren 

1 0 S ignal verfalschungen. 

Sicher vermeiden lassen sich die Fehler dadurch, daQ die Lange L der Filterimpulsancwort, die den 
Bewertungsfaktoren entspricht, geeignet begrenzt wird. Nachdem aber bei sehr vielen adaptiven 
Merkmalsanderungsverfahren diese Bewertungsfaktoren nur im Frequenzbereich geliefert und fort- 
laufend angepaBt werden, ware nach jeder Anderung ein Filter-Entwurfsverfahren, z.B. nach der Stra- 

15 tegie des "Frequency-Sampling" (siehe Z.B. das o.g. Buch von Oppenheim und Schafer), durchzu- 
fuhren. Dieses erforderte viel Rechenzeit, wurde somit die fortlaufende Verarbeicung unterbrechen 
und die Echtzeitanwendung der meisten der o.g. Merkmalsanderungsverfahren unmoglich machen. 

Die Erfindung lost diese Probleme dadurch, daQ in das adaptive "Overlap-Add" -Verfahren der im 
20 kennzeichnenden Teil des Anspruchs I beschriebene Zusatzalgorithmus eingebracht wird, der es er- 
laubt, die "Time -domain-aliasing" -Fehler sicher unter einem wahlbaren Grenzwert zu halten. Der 
Rechenaufwand fur diesen Zusatzalgorithmus ist gering, er macht in den meisten Fallen nur einen 
kleinen Bruchteil des fur das "Overlap-Add" -Verfahren ohnehin notigen Aufwandes aus. Des weiteren 
hat es sich als vorteilhaft erwiesen, daQ der Rechenaufwand des erflndungsgemaCen Zusatzalgo- 
25 rithmus 1 abnimmt, wenn groQere Fehler zugelassen werden. Dadurch kann fur jede Anwendung ein 
optimaler KompromiB zwischen Rechenaufwand und Qualitat des Ausgangssignals gefunden werden. 
Das Konzept und die Anwendung des erfindungsgemaQen Feblerbegrenzungsverfahrens werden nach- 
stehend anhand der Figuren beispielhaft erlautert, ohne dafl dadurch andere, fur fachkundige Personen 
naheliegende Ausfuhrungsformen ausgeschlossen oder eingeschrankt werden sollen. 

30 

Es zetgen 

Fig. 1 : das Blockschaltbild eines Verfahrens zur Merkmalsanderung gemafl dem herkommlichen, 
im Oberbegriff des Anspruchs I beschriebenen Stand der Technik. 

Fig. 2: ein beispielhaftes Schema zur Rekonstruktion des Ausgangssignals durch uberlappende 
Addition von Signalabschnitten. 

Fig. 3: das Blockschaltbild eines erfindungsgemaOen Verfahrens zur Merkmalsanderung. 

Fig. 4: die erfindungsgemaOen Schieberegister 131 und 134 mit eingetragenen beispieihaften 
Frequenzgang- und Fensterfolgen. 
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* 5: das Ergebnis der RuckrransfonBation der in der Fig.4 dargesreilten Fensterfixnktioa ir 
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D» F* 1 ze.gr das Biocfccnaltbild eines Cerates zur Merkmalsandenzng, welches dem i m Obe, 
begnff des Ansprucbs 1 beschriebenen Stand der Technik entspricnt. Dabei wird angenommen. da 0 
das Emgangssignal x(t), welches Qber die Leitung 10 an den Analog-Digital-Wandier 20 herangefuhn 
w«i, emen kontinuierlicaen elektriscbeu Spannungsverlauf uber der Zeit t darstellt. In bestitnmten 
Anwendungsrallen wurde dieses Signal durch einen geeigneten Wandler (z.B. Mikrofon) aus dem 
Verfauf emer anderen physikalischen GroBe (z.B. Schalldruck) in bekannter Weise erzeugt. Der Block 
20 (Aaaiog-Digical-Wandler) beinhaltec auch die Abcast-Halte-Funkcion, wenn eine solche fur das 
Wandlungsverfahren erforderiich ist. Dabei wird angenommen, daO die Wene mit gentigender Ge- 
nauxgkezt, dh. genagender Anzafal Bits je Wert x(n), quantisierc werden. Ober die Leitung 30 wird 
-nut one 2e itliche Foige x(n) digitalisierter Abtastwerre des Eingangssignals an das Scfaieberegister 
40 und an die Verarbeitungseinneit 50, weiche die Adaptationsstrategie realisien, herangefuhrt. Dabei 
stellt n dea Zihlindex der aufeinanderfolgenden Abtastintervalle dar. Die Cbertragung kann wahl- 
wetse parallel oder bitseriell erfolgen; dies ist fur das hier zu erlautemde Verarbeitungsverfahren on- 
ernebhch. Im Scfaieberegister 40 werden M Abtastwene x(n) aufbewann. Wird ein neuer Wen ein*e- 
lesen, so wird der aiteste gespeichene Wen binausgeschoben und geht verlorer, Sind K neue Wene 
euxgelesen, beginnt ein Verarbeitungszyklus. M mufl ein ganzzahliges Vielfaches von K sein, M=K ist 
moghcb, aber unzweckmaflig. Dem Scfaieberegister 40 werden M Abtastwene entnommen und uber 
d:e Lextung 60 einer Multipliziereranordnung 70 zugefuhn. Die Leitung 60 ist als Bandleiron* aus M 
parallelen Emzelieitungen dargestellt, um zu symbolisiereu, dafl in jeweils einem Verarbeitungsschritt 
em Block von M Abtastwenen bearbeitet wird. In der Multipiziereranordnung 70 wird jeder Ein- 
gangswen mit je einem Wen einer Fensterfunktion w(n) multiplizien, die in einem Speicher 80 
berextgehalten wird. Der Block 70 kann M parallel arbeitende Multiplier enthalten. Es kann aber 
auch em Multzpiizierer geuugen, wenn er schnell genug ist, urn innerhalb von K Abtastintervallen M 
Multzplzkationen in Zeitmultiplex-Technik abzuarbeiten. Die Fensterfunktion w(n) mufl so gewahlt 
werden. daB sie bei Verschiebung urn Vielfache von K und Aufaddition eine konstante Folge liefen 
wie , n der genannten Veroffentiichung von Allen und Rabiner naher erlauten wird 
D.e M fensterbeweneten Eingangsdaten werden erganzt durch N-M Nullwene. was durch die Anord- 
nung 90 symbolisiert wird. Der so entstandene Datenblock aus M Wenen wird einer diskreten Spek- 
traltransformation 100 zugefuhn. Allen und Rabiner haben hierfur die diskrete Fouriertransformation, 
« der recheneffizienten Ausfuhrungsform der "Fast Fourier Transform" (FFT) vorgeschlagen. Alter- 
nate Realisierungsmaglichkeiten fur die Spektraltransformation 100 sind in den Anspnichen 5 bis 8 
angegeoen. Die Spektraltransformation liefer, ein Spektrum X(v) aus N diskreten, in vieien Falle- 
komplexen Wenen, welcfaes uber die Leitung 110 einer Multipliziereranordnung 120 zugefuhn wird 
Don erfoigt die Merkmalsanderung, indem jeder Wen X(v) mit einem Wen der im Speicher 130 

bereugehaltenen Frequenzgangfunktion H(v) multiplizien wird. Es ergibt sich ein diskretes Aus«an*s- 
spekrrum 3 3 
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Y(y) = X(y).H(y) fcr v =0,1,2,.. .JV-l. (l) 
Die durch den Block 50 symbolisierte Adaptatiatisstrategie andert von Zeit zu Zeit einzelne oder alle 
Werre von H(n) im Speicher 130 uber die Leitung 140, abhangig von Verandenmgen im Eingangs- 
signal x(n). Weil die Adaptationsstracegie niche Gegenstand der vorliegenden Erfindung ist, wird diese 
hier niche naher beschrieben. 
5 Das Ausgangsspektrum Y(n) wird uber die Leicung 150 einer inversen Spektraltransfonnation 160 
zugefiihrt, welche die Umkehrung der im Block 100 verwendeten Transformation darstellt, so daO N 
digitale Abtastwerte yj(n) eines Zeitbereichssignals errechnet werden. Dabei bezeichnet i den laufen- 
den Index der Verarbeitungsschritte. Wie oben beschrieben, werden in jedem Verarbeicungsschritt aus 
dem Schieberegister 40 nur K neue Abtastwerte, aber M-K altere Abtastwerte, welche bereits einmal 

10 oder sogar mehnnals in fniheren Verarbeitungsscbritten verwendet wurden, der Verarbeicung zuge- 
fiihrt. Die Verarbeicung erfolgt also in uberlappenden Datenabschnitten. Dementsprechend ergibt sich 
am Ausgang die endgultige Ergebnisfunktion y(n) durch zeitverschobenes Aufaddieren mehrerer 
Signalabschnitte yi(n). Dies wird durch die Addiereranordnung 170, die beiden Schieberegister 180 
und 190 sowie den Umschalter 200 bewirkt. Die Addiereranordnung 170 fuhrt in jedem Verarbei- 

15 tungsschritl N-K Additionen aus. Dabei ist es gleichgaltig, ob N-IC paralleie Addierer realisiert sind 
oder ob ein schneller Addierer die Additionen im Zeitmultiplexverfahren abarbeitet. Die Ergebnisse 
der Additionen und die K neuesten, noch nicht der Addition unterworfenen Ergebniswerce werden 
parallel ins Schieberegister 180 eingespeichert und dann seriell nach oben herausgeschoben. Dabei ist 
der Umschalter 200 zunachst nach rechts geiegt, so dafl die K altesten im Schieberegister 180 gespei- 

20 cherten Werte dem Digital-Analog-Wandler 210 zugefuhrc werden und somit einen Teil der Aus- 
gangsfolge y(n) bilden. Dann wird der Umschalter 200 nach links umgeschaltet und die restlichen N-K 
Werte werden ins Schieberegister 190 umgespeichert, so daD diese fur weitere Additionen zur Ver- 
fugung stehen. 

Die Fig. 2 zeigt fQr die beispielhaften Werte N=256 und K=64 die zeitverschobene und uberlappende 
25 Aufaddition der Signalabschnitte y^n), i=0,l,2,... t zur Ausgangsfolge y(n). Dabei sind untereinander- 
stehende Werte der Teilfolgen yj(n) zu addieren. 

Der in der Fig- 1 dargestellte Digital-Analog-Wandler 210 erzeugt aus y(n) ein kontinuierliches elek- 
trisches Ausgangssignal y(t). Dieses kann uber die Leitung 220 einer Weiterbearbeitung zugefiihrt 
werden, beispielsweise auch einer Umwandlung in eine andere physikaiische Darstellungsform. 

30 

Die in der Fig. 1 dargestellte Verarbeicung erfolgt aber nur dann frei von M Time-domain-aiiasing"- 
Fehlern, wenn zu der im Speicher 130 vorhandenen Frequenzgangfunktion H(v) im Zeitbereich eine 
Impulsantwort h(n) gehort, welche nicht mehr als L von Null verschiedene Werte aufweist und wenn 
L die Bedingung 

L<N-M-»-l (2) 
35 erfullt. Wenn die Adaptationsstrategie 50 fortlaufend neue Frequenzgangfiinktionen H(v) bereiKteilt, 
ist die Uberprufung der Bedingung (2) so aufwendig, daB Echtzeitverarbeitung im allgemeinen nicht 
mehr moglich ist, wie oben bereits dargeiegt wurde. 
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Schieberegisrerausgangen ansteheaden Werte H(-B) und G(-B) werden vom Multiplizierer 133 
multipliziert, der Addierer 136 addiert Null, so daG das Produkt unverandert in den Speicher 137 
gelangt. Nun werden beide Schieberegister, 131 und 134, einmal geschoben. Der Multiplizierer 133 
bildet das Produkt H{-B + 1) • G(-B + 1) . Dieses wird im Addierer 136 zum gespeicherten Produkt 
5 H(-B)-G(-B) addiert und die Summe im Speicher 137 aufbewahrt. In dieser Weise wird 
fortgefahren, bis J Teilprodukte H(y)-G(y) addiert und die Gesamtsumme im Speicher 137 
eingespeichert ist. Nun steht am Ausgang des Speichers 137 der erste Wert des durch die Faltung 
modifizierten Frequenzgangs H(0) an. Durch einmaliges Schieben des hier als Schieberegiscer 
ausgebildeten Speichers 130 wird dieser Wert im Schieberegister 130 eingespeichert. Multiplizierer 

10 133 und Schieberegister 134 werden daraufhin angehalten, wahrend das Schieberegister 131 N-J+L 
weitere Schiebeschritte ausfiihrt. Damit hat letzteres Register seit Beginn des Falrungsvorgangs N+l 
Schiebeschritte durchgefuhrt, wodurch sich eine urn einen Wert kieinere Verschiebung von H(v) als in 
Fig. 4 dargestellt ergibt, d.h. H(-B+l), hier H(-3), steht nun in der Ausgangszelle des Schieberegisters 
131. Wenn nachfolgend der Speicher 137 -geioscht wird und anschlieOend das Produkt 

15 H(-B + 1) • G(-2?) in den Speicher 137 eingespeichert wird, so ist dies der erste Schritt zur 
Berechnung des zweiten modifizierten Frequenzgangwens H (I) . Wenn die Berechnung dieses 
Wertes vollendet ist, wird auch dieser Wert aus dem Speicher 137 in das Schieberegister 130 
geschoben. In dieser Weise wird fortgefahren, bis schlieQlich N modifizierte Frequenzgangwerte H(v) 
im Schieberegister 130 stehen. Von dort konnen diese wiederhoit parallel in die Multipiizierer- 

20 anordnung 120 ubertragen werden, urn die Merkmalsanderung bei den Eingangsdatenblocken so zu 
bewirken, wie das oben bereits beschrieben wurde. Erst dann, wenn die Adaptationsstrategie 50 fest- 
legt, dafl eine neue Frequenzgangfunktion H(v) zur Anwendung kommen muB, wird die erfindungs- 
gemaCe Fensterung wieder aktiviert. 

25 Bei der beispieihaften Anwendung eines gemafl der Fig. 3 realisierten Gerates zur Dynamikkompres- 
sion von Sprachsignalen zeigte es sich, daB ein solches Gerat kiein, leicht und daher als experimentelle 
Horhiife gut tragbar realisien werden kann, wenn die Parameter wie folgt spezifiziert werden: Ein- 
gangssignal-Abschnitt M=180; Anzahl neuer Werte in jedem Eingangssignal-Abschnitt: K=90; Aus- 
bildung der Transformation als FFT der Lange N=256; Lange des Fensters G(v): J=9. Bei der Reali- 

30 sierung mit Hiife des Assembler-Codes eines modernen schnellen Signalprozessors, z.B. des 
DSP56L002 der Motorola Semiconductor Ltd., erwies sich der Rechenaufwand des erfindungs- 
gemaDen Gerates (gemafl Fig. 3) gegentiber einem nach dem konventionelien Stand der Technik her- 
gestellten und der Fig. I entsprechenden Geriit als nur geringfugig, namlich urn ca. 12% erhoht. 

* 35 Aus theoretischen Abhandlungen, Z.B. in 'The digital prolate spheroidal window" von T. Verma, 3. 

Bilbao und T.H.Y. Meng, erschienen in den Proceedings of the International Conference on Acoustics, 
Speech and Signal Processing (ICASSP) 1996, veranstaltet vom IEEE in Atlanta/USA, Seiten 1351- 
1354, konnte gefolgert werden, daD die geringsten verbleibenden "Time-domain-aliasing"-Fehler mit 
der "Prolate Spheroidal" Fensterfunktion erzielt werden konnten. In der Fig. 4 ist im Schieberegister- 
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Patentanspriiche 

1. Verfahren zur Merkmalsanderung bei eindimensionalen, in eine digitaie Darsteilung 
umgewandelten Signaien unter Verwendung des adaptiven "Overlap-Ad^'-Algorithmus', welcher 
die Merkmalsanderung nach einer geeigneten diskreten Spekcraltransformarion durch Multiplica- 
tion im Frequenzbereich durchfuhrt und anschlieflend das Ausgangssignal durch entsprecbende 
inverse diskrete Spektraltransformation sowie durch uberlappende und verschobene Addition 
mehrerer, von der inversen Spektraltransformation gelieferter Signaiabschnitte erzeugt, dadurch 
gekennzeichnet, dafl vor den Multiplikationen im Frequenzbereich die Frequenzgangfunktion mit 
einer wahlbaren diskreten Fensterfunktion, welche erheblich geringere Lange aufweist als die 
Frequenzgangfunktion, gefaltet wird. 

2. Verfahren nach Anspruch 1, dadurch gekennzeichnet, dafl die Fensterfunktion im Frequenzbereich 
nur positive Werte und einen glatten Verlauf aufweist, im Zeitbereich aber eine ausgepragte 
Strukcur mit einer groflen positiven Hauptkeule und alternierenden, dem Betrage nach kleinen 
Nebenmaxima und Nebenminima besitzt. 

3. Verfahren nach Anspruch 1, dadurch gekennzeichnet, dafl die Spektraltransformation als diskrete 
Fouriertransformation und die inverse Spektraltransformation als inverse diskrete Fouriertrans- 
formation ausbildbar ist. 

4. Verfahren nach Anspruch 3, dadurch gekennzeichnet, dafl die diskrete Fouriertransformation und 
die inverse diskrete Fouriertransformation mit Hilfe des Algorithmus' der "Fast Fourier 
Transform" (FFT) durchgeftihrt werden. 

5. Verfahren nach Anspruch 1, dadurch gekennzeichnet, dafl die Spektraltransformation als diskrete 
Cosinus-Transformation ("Cosine Transform") und die inverse Spektraltransformation als inverse 
diskrete Cosinus-Transformation ausgebiidet ist. 

6. Verfahren nach Anspruch 1, dadurch gekennzeichnet, dafl die Spektraltransformation als diskrete 
Haar-Trans formation und die inverse Spektraltransformation als inverse diskrete Haar- 
Transformation ausgebiidet ist. 

7. Verfahren nach Anspruch 1, dadurch gekennzeichnet, dafl die Spektraltransformation als diskrete 
Walsh-Hadamard-Transformation und die inverse Spektraltransformation als inverse diskrete 
Walsh-Hadamard-Transformation ausgebiidet ist. 

8. Verfahren nach Anspruch I, dadurch gekennzeichnet, dafl die Spektraltransformation als diskrete 
Hartley-Transformation und die inverse Spektraltransformation als inverse diskrete Hartley- 
Transformation ausgebiidet ist. 
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18. Gerat zur Durchfuhrung des Verfahrens nach einem oder mefareren der Anspruche 1 bis 13 , 
dadurch gekennzeichnet, daQ das Verfahreu als eine Folge von Assemblerbefehlen dargestellt 
wird, welche auf einem Mikrocomputer oder einem spezialisierten Signalprozessor ausfiihrbar ist. 

19. Gerat nach Anspruch 18, dadurch gekennzeichnet, dafl die Folge von Assemblerbefehlen durch 
ein Compilerprogramm erzeugt wird, welches die dazu benotigte Information einem in einer 
hoheren Programmiersprache geschriebenen Programm enmimmt. 
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